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Abstract. The basic properties of traffic flow are analyzed using a sim- 
ple deterministic one dimensional car following model with continuous 
variables based on a model introduced by Nagel and Herrmann [Physica 
A 199 254-269 (1993)] including a few modifications. As a first case we 
investigate the creation and propagation of jams in a platoon generated 
by a slow leading vehicle. In a second case we look at a system with the 
size L, periodic boundary conditions and identical vehicles. A strong de- 
pendence on the initial configuration of the fundamental diagram's shape 
can be found. 



1 Definition of the Model 

To get a parallel update of all vehicles, the first step is to change the velocity of 
all vehicles taking into account the threshold parameters a and [3 according to 
the rules defined below. In the second step the position of all vehicles is changed 
using the velocity calculated in the first step. 

Position, velocity and threshold parameters of the vehicles are continuous 
variables. At is the time step used in the simulations. 

1st step - velocity: A vehicle decelerates if the headway distance Axi{t) is 
smaller than a safety distance a. The headway distance after the deceleration 
step is determined by 5. 

Axi(t) - Vi{t)At < a v,{t ^ At) ^ max ^0, — (1) 

The acceleration of a vehicle depends on a threshold parameter (3 and the max- 
imum velocity v„iax- 

Ax,{t) - v,{t)At > (3^ v,{t + At) = min {v.max, v^) + aAt) (2) 

The acceleration coefficient is determined by the headway distance if 7 is greater 
than /3 and the headway distance is smaller than 7: 

/ Ax,{t) \ 

a = Umax max 1, . (3) 
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2nd step - position: The positions are changed. 

x,{t + At) ^ x^{t) +v{t + At)At (4) 

Using a constant acceleration coefficient (7 < /?) we can compare this model 
to the deterministic case of the Nagel-Schreckenberg model [|| , as the velocity is 
determined by de- and acceleration with the acceleration causing a discretisation 
of the velocity. 

2 Platoon 
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Fig. 1. The figures on the left side show the trajectories of vehicles and jams 
of diff'erent sizes. On the right side the oscillations of the velocity and relative 
position (position in the coordinate system of the leading vehicle) of the second 
and third vehicle are shown including the threshold conditions for accelerating 
and breaking. Parameters: a = 15,/3 = 35,7 = 100,(5 = 20, Umax = ^,Vmax = 

30,Vlead = 25 

A slow leading vehicle that is followed by faster ones creates a platoon of 
vehicles. In this platoon, jams are generated directly behind the leading vehicle. 
They propagate either through the entire platoon or at least through parts of it. 
These jams are created because a vehicle is not able to slow down exactly to the 



3 



velocity of the one ahead. In our model this is forced by the parameter S, that is 
the headway distance after a breaking step. For 5 > a the headway distance a at 
the beginning of the breaking step is smaller than the headway distance 6 after 
the breaking. The consequence is, that the vehicle has to increase the headway 
distance within the breaking time step by decelerating to a slower velocity than 
the car ahead. This leads to an oscillation in the velocities and the headway 
distances of the vehicles and can add up in a constructive way generating jams 
with the size of the entire platoon. 

3 Periodic Boundeiry Conditions 

In our study, we work on the fundamental diagram (flux vs. density) of a sys- 
tem with the size L and periodic boundary conditions using identical vehicles. 
Strong dependencies on the fluctuations in the initial configuration have been 
found and investigations of three different cases have been made. On one side we 
look at a highly symmetric case without any fluctuations using equidistant start 
positions. On the other side we use homogeneous random start positions. Be- 
tween these two cases, random fluctuations are used to move the vehicles out of 
their equidistant positions. The maximum of these fluctuations, however, is con- 
trolled by a parameter AL that can be in the range from zero up to the system 
size. The idea is that with the parameter AL the fluctuations or the disorder of 
the initial state can be characterized including the two extreme cases mentioned 
above representing the configurations for the minimum and maximum values of 
AL. For these extreme cases analytical solutions can also be found. The velocity 
of all vehicles is set to zero in the initial configuration and the same value is used 
for the parameters a and S. 

3.1 Equidistant Starting Positions 
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Fig. 2. Flux for different simulation time-steps At and equidistant start posi- 
tions. Parameters: a = ^ = 15, /? = 35, 7 = 10, Umax = 1, Vmax = 30 
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In this highly ordered case the model generates no jams and the fundamental 
diagram is a straight line with a slope of Vmax- The maximum flux can be 
obtained from the accelerating condition of the model. 

Pcrit = -TT ; ^ (5) 



^{Pcrit) = Vmax Pent (6) 

For At — > the maximum flux jumps to zero at the maximum density Pmax = 
1//3, which is determined only by the acceleration threshold parameter /J. 



3.2 Homogeneous Random Starting Positions 




Fig. 3. The left side shows density waves in a space time diagram. The right 
side shows the flux for random initial positions and various safety distances. 
Parameters: a = (5 = 15, /3 = 35, 7 = 10, Umax = 1, Vmax = 30 



In this case we flnd a triangular-shaped fundamental diagram with a maxi- 
mum density Pmax = l/a and a critical density Pcrit with maximum flux ^max- 
Up to the critical density the system is free-flowing and from the critical density 
up to the maximum density there are both, jammed and free flowing areas. Be- 
cause of the sharp boundaries between congested and free areas, the system can 
be seen as a mixture of two phases like water and ice. The size of the jammed and 
the free phase, respectively, can be obtained through the normalized (divided by 
the maximum flux) falling straight line of the fundamental diagram. 

Looking at the trajectories of vehicles entering a congested area, it can be 
proven that the density of the congested area depends on the velocity of the entire 
jammed area itself. In the final state of a finite system, only the congested areas 
with the slowest velocity survive. Taking this into account and neglecting the 
acceleration and deceleration stripes, there are two phases characterized by the 
density and the velocity of the vehicles within each phase. In order to obtain the 
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ftindamental diagram for this two-phase system the mean vefocity is calculated. 

^ - ^ E - ^ - - ^ (^) + P) 

The indices j and / are used for jammed and free phase, respectively. 

The density of the congested area can be obtained from the braking rule of 
the model. 

Pj = ^—ir (8) 

The outflow of a congested area determines the density of the free flowing area 
that is equal to the critical density. As a result of analyzing the trajectories, it 
is found that the time Z\T, that is needed to accelerate and pass the window 
/3 — a, is an important value in this model. 



AT 



1 /I 2(/3- 



2 V 4 aZ\<2 



1 At (9) 



The braces [..J operate on the given argument by rounding down to the next 
integer value. The result for the critical density using v^aax = is therefore: 

1 , ^ 

Pent ^ Pi = \ XT— 7 TTTTT (10) 

a^VjAt+{Vraax-Vj)AT 

From equation (|^, ||, H, |l^) the exact solution for the flux using a constant 
acceleration coefficient can be obtained: 

m - ^{p)P = ^+p (^'. - ^^^) (11) 

Investigations on the connection between the jam velocity and the slope of the 
fundamental diagram have shown an exact agreement with the results of Lighthill 
and Whitham js). 

3.3 Equidistant Starting Positions with Random Fluctuation 

Equidistant positions including small fluctuations are used for the initial config- 
uration. 

xdt = 0) = - + ^ (13) 
P 2p 

The random number p generates the fluctuations and characterizes their magni- 
tude with respect to the density. 
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Fig. 4. The left side shows the trajectories of vehicles and jams with a velocity 
greater than and a density of 20/km. The fundamental diagram is shown on 
the right side. Parameters: a = ^ = 15, /3 = 35, 7 = 10, amax = 1, t^max = 30,p = 
0.75 



A higher maximum flux and critical density are obtained compared to the 
homogeneous random case. Prom the critical density up to the maximum equidis- 
tant density wc get congested areas consisting of vehicles with a velocity 
Vi> Q that leads to a greater mean velocity and a greater flux. From the equidis- 
tant maximum density 1//? up to a cutoff density p^ut, the flux matches exactly 
the homogeneous random case. The cutoff density pcut can be calculated directly 
from the acceleration condition including the fluctuations. 

Pent = ^ (14) 



4 Outlook 

A simple Id traffic model with deterministic propagation rules has been studied. 
A strong dependence on the fluctuations of the initial state has been found and 
for the extreme cases analytical solutions can be obtained for the fundamental 
diagrams. 

It would be interesting to generalize the present model to more than one 
lane and to introduce more realistic velocity dependent threshold parameters. In 
the vicinity of the critical density investigations with respect to bistable states 
would be highly interesting. 
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